Nonequilibrium fluctuations of an interface under shear 
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The steady state properties of an interface in a stationary Couette flow are addressed within the 
framework of fluctuating hydrodynamics. Our study reveals that thermal fluctuations are driven 
out of equilibrium by an effective shear rate that differs from the applied one. In agreement with 
experiments, we find that the mean square displacement of the interface is strongly reduced by the 
flow. We also show that nonequilibrium fluctuations present a certain degree of universality in the 
sense that all features of the fluids can be factorized into a single control parameter. Finally, the 
results are discussed in the light of recent experimental and numerical studies. 
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I. INTRODUCTION 

Soft matter systems driven far from equilibrium by a 
shear flow manifest striking properties as a resuh of their 
sensitivity to external fields For instance, hydrody- 
namics may either enhance or suppress coarsening pro- 
cesses i^. Coupling with the flow can also lead to the 
emergence of new, shear-induced phases that are exclu- 
sively out-of-equilibrium structures 0, whereas spa- 
tiotemporal oscillations and rheochaos are observed in 
shear-banding systems d, Q . 

Despite substantial progress, a fundamental under- 
standing of complex fluids under shear remains a chal- 
lenging question. The flrst reason lies in the nature of 
the coupling between structure and flow. Since the lo- 
cal structure of soft materials is readily reorganized by 
an external flow, it has in turn a signiflcant impact on 
the flow itself. The complexity of this feedback mecha- 
nism makes that most theoretical studies are based - to 
a variable extent - on phenomenological models fH, 
The second point is that our understanding of nonequi- 
librium steady states (NESS) is yet far from complete. 
Recently, there have been several attempts to construct 
a unified framework to describe equilibrium and nonequi- 
librium phenomena Q. In particular, fluctuations theo- 
rems [8] or extended fluctuation-dissipation relations [9J 
have been suggested for complex fluids under shear. But 
these relations still remain at the conceptual level and 
the link with experimental observables - power spectra, 
correlation functions - has not been clarified yet. 

In this paper, we investigate the statistical properties 
of a liquid-liquid interface in a stationary flow. We ar- 
gue that this system is sophisticated enough to capture 
the relevant features of NESS but remains simple enough 
to be addressed analytically. Applications of this issue 
might be expected for instance in microfluidics, since in- 
terfacial phenomena are inevitably enhanced as the size 
of the system is reduced. Still, the question of interface 
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fluctuations under shear has received only little atten- 
tion so far. From a fundamental viewpoint, equilibrium 
properties of liquid-liquid interfaces are now well estab- 
lished [To| , and experimental as well as numerical studies 
have validated the capillary wave model down to almost 
molecular scales [HI [12] . This topic has experienced a 
renewed interest in recent years with the discovery that 
phase-separated colloid-polymer mixtures may have an 
extremely low interfacial tension [13^ . Interface fluctua- 
tions can then be analyzed in real time and space using 
elementary video-microscopy techniques. The versatility 
of this method has allowed Derks and collaborators to 
study the statistical properties of an interface exposed 
to a shear flow fl3|. They found that the coupling with 
the flow leads to a strong reduction of thermal fluctu- 
ations, while the correlation length increases. But this 
second point is in disagreement with recent Monte Carlo 
simulations of a driven Ising model , therefore raising 
the fundamental question of what features of interfaces 
under shear are actually universal. 

Statistical fluctuations of an interface are driven by the 
random forces that spontaneously occur in the bulk [l6j . 
To describe nonequilibrium properties, the main issue is 
thus to properly account for the coupling between the 
bulk and the interface. This can be achieved on the basis 
of fluctuating hydrodynamics (FH) Indeed, FH has 
been successfully applied to various NESS situations, and 
experimental validations have been obtained, e.g., for flu- 
ids in a temperature gradient [3. The purpose of this 
paper is to apply FH to an interface under shear. 

We shall proceed as follow. In Sec. |TT1 we explain how 
an equation of motion for the interface can be obtained 
from FH. We show that the distortion of capillary waves 
by the flow leads to a mode-coupling equation that is 
discussed in Sec. IIIII This allows us to extract NESS 
properties under a stationary flow in Sec. IIVI In partic- 
ular, we find that the fiuctuations are smoothed out by 
the flow. The results are then discussed in Sec|V] in the 
light of recent experimental and numerical data available 
in the literature. Finally, we conclude the paper with 
a short summary of our results. For the sake of clarity, 
details of the algebra are presented in Apps. [XHdI 
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II. HYDRODYNAMIC FORMULATION 

We first set up a hydrodynamic theory to account for 
the coupUng between the surface and the bulk. Following 
the usual hypothesis of capillary wave theory, we assume 
the existence of an intrinsic interface separating two im- 
miscible fluids. For moderate deformations around the 
xOy horizontal plane, the position of the interface can 
be described by a single valued function z — h{x,y,t). 
Along this article, properties of the upper (resp. lower) 
fluid are labeled with the subscript i — 1 (resp. i — 2). 
Each phase is characterized by its mass density pi and its 
viscosity rji. We also define fj = (r/i -|-ry2)/2 the mean vis- 
cosity and Ap = p2 — pi > the mass density difference. 
The surface is further characterized by the interfacial ten- 
sion a and the capillary length 1^ = y^a/{Apg), with g 
the gravitational acceleration. 

The system is schematically drawn in Fig. [T] The av- 
erage position of the interface is z = 0. It is confined 
between to walls, the thickness of each fluid layer be- 
ing Li and L2 with L = Li + L2. A planar Couette 
flow is induced by moving the walls of the shear cell at 
constant velocity along the x direction. We define x, y 
and z respectively as the velocity, vorticity and velocity 
gradient directions. Assuming that the no-slip condition 
applies on the walls of the cell, the fluid velocity v sat- 
isfies Y{x,y,Li) = ViBx and v(a;,?/,— L2) = —V2^x- The 
total shear rate is then 

where we define 7^ = Vi/Li the shear rate in each phase. 
Without loss of generality, we assume in the following 
that Vi , Li , V2 and L2 are chosen so that the plane of zero 
shear coincides with the average position of the interface. 
Other situations can be deduced thanks to a Galilean 
transformation. 

For usual fluids at room temperature, thermal fluctua- 
tions occur in the overdamped regime of capillary waves. 
In order to describe nonequilibrium effects, the analysis 
is performed within the framework of fluctuating hydro- 
dynamics [13, The starting point is the stochastic 
version of the Stokes equation 

77,V2v - Vp + Pig + V • s = , (2) 

with V the velocity field, p the pressure, and g = —gez- 
Eq. ^ is solved together with the incompressibility con- 
dition 

V-v = 0. (3) 

Thermal fluctuations are accounted for through the ran- 
dom part of the stress tensor s. Its components s^i, 
(with fj.jh' = X, y, or z) are stochastic forces that stem 
from the microscopic degrees of freedom of the fluids. 
Close to equilibrium their correlations are given by the 
fluctuation-dissipation theorem, but such a relation is not 
expected to hold beyond the regime of linear response. 
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FIG. 1: Schematic representation of the system. The unit 
vector n is normal to the interface and is pointing towards 
the upper fluid. 

Here however, we shall take advantage of the separation 
of time scales between the collective modes under study 
- the fluctuations of the interface - and the molecular 
scales of the heat bath - the fluid constituants. The re- 
laxation of an interface is characterized by the capillary 
time Tc — 2fjlc/(T; it ranges from milliseconds for usual 
interfaces (a « 10~^ mN/m) to a few seconds for ultra- 
soft interfaces (tr w 10~^ mN/m). On the other hand, 
a particle of fluid on either side of the interface diffuses 
over its own diameter a on a time-scale Tb = rjiU^ /{kBT). 
This time scale is at most of the order of Tf, w 10"'^ s for 
particulate fluids with a « 100 nm [3], but it is several 
orders of magnitude smaller for molecular fluids. Since 
the relevant regime considered in this work corresponds 
to jTc 1, the applied shear rate is too small to signif- 
icantly affect the thermal motion of individual particles. 
In this small Peclet number limit {jri, ^1), the statis- 
tics of the heat bath is not affected by the flow and we 
can assume that the stochastic variables have zero mean 
value and correlations given by [l3, [3 

(s^^ {r,t) Sf_,'^' (r',i')) ^ '^ksTiji {S^^'S^^^ +S^^'S^^') 

xS{r-r')6it-t') ,(4) 

with fcs the Boltzmann's constant and T the tempera- 
ture. 

Eqs. ([2]) and ([3]) are then solved above and below the 
interface, and the solutions matched with the appropriate 
boundary conditions. The latter have to be enforced at 
the instantaneous location of the interface z = h{x,y,t). 
Explicitly, we need to express the continuity of the ve- 
locity 

M, = o, (5) 

as well as the continuity of the stress 

[T]^.n = an(V-n) , (6) 
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with the notation [f]^^^ = /(^o~) — /{zq ), the hmit be- 
ing taken respectively from above and from below. In 
Eq. (O, T = t + s is the total stress tensor [H, [l^. The 
components of t read = ^P^^i.v + {d^Vv + d^v^), 
with fijV — X, y, or z; the random part s is defined in 
Eq. ([2]). The unit vector n is normal to the surface, 
pointing towards the upper fluid. It depends on the local 
conformation of the interface 1201 
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VT+TW\ 1 



-dxh\ 
dyh 



(7) 



Finally, once the velocity field is fully characterized, an 
equation of motion is obtained thanks to the kinematic 
relation 



- V|| • V||/i = 



(8) 



the velocity being evaluated at z — h. In this equation, 
the parallel components of a vector field f — {f^, fy, fz) 
are defined as f|| — [fx, fy)- For the sake of completeness, 
the derivation of the boundary conditions ^ and ([8]) is 
reminded in Appendix El 



III. THE INTERFACE EQUATION 

It has been appreciated for a long time that intcrfacial 
fluctuations can be distorted by an external flow 21,, J2J . 
The description of this effect requires to go beyond the 
linear analysis generally used to derive the dispersion re- 
lation for capillary waves. To proceed, we follow a re- 
cursive scheme that has proved its worth, e.g ., in the 
context of polymer-membrane interactions [23|. We as- 
sume that the deformation can be written h{x,y,t) = 
eu(x,y,t), with u(x,y,t) ~ 0(1)- The dimensionless 
parameter that governs this small-gradient expansion is 
£ — yJkBT/{<7l1). The velocity and pressure fields are 
then expressed as 
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The lowest-order term is simply the Couette flow solution 
for a planar interface: v('''(r) = 712:62; if z > 0, and 
v(o)(r) =72ze^ if z < 0. 

For n > 1, the linearity of the Stokes equation implies 
that each term of the series obeys Eq. ^ . The coupling 
between successive orders originates from the boundary 
conditions. Indeed, the latter have to be enforced at the 
position of the interface. It is therefore expected that the 
Taylor expansion of the n*''-order field, when evaluated at 
z = eu{x, y, t), involves contributions of all orders k > n. 
To be more explicit, consider for instance the velocity 
field. Up to second order, it is given by 

v{eu) = vW(0)+e [v(i)(0) + ua^v(°)(0) 



v(2) (0) -f u9,v(i) (0) + y a,V°) (0) 



[Note that we have only written the z-dependence of 
■v{x,y, z,t)]. We proceed likewise for all fields in Eq. ([5]) 
- ([8]). Clearly, the Taylor expansion of the boundary con- 
ditions involves several second-order contributions, mak- 
ing the algebra quite cumbersome. Moreover, the study 
of the fluctuations requires us to solve the full three- 
dimensional problem. We thus defer the technical details 
to App.lBlandlCl and we now focus the discussion on the 
main outcomes. 

Since the problem is invariant by translation parallel 
to the horizontal plane, it is natural to switch to the (2D) 
Fourier representation /i(q, t) = / d^rpexp[— zq-r]ft,(r, f), 
with r = {x, y) and q = {qx,qy)- We also define the norm 
of the wave vector q — |q| = {q^ + q^Y^'^ ■ After some 
algebra, we find that the relaxation of a fiuctuation mode 
with wave vector q follows a mode- coupling equation 

dth = -;^/i(q, t) - i'^eff J ^;^kxh(k, i)/i(q - k, t) 

+ ip{<l,t). (9) 

This equation constitutes the first main result of this pa- 
per. It involves a number of contributions that we now 
discuss. First, Tq = Af]q/{(j{q^ + is the equihbrium 

relaxation time of the interface. This result shows that 
there is no direct coupling at linear order (at least in the 
viscous regime). Advection of the interface by the flow 
takes the form of a convolution between all modes with 
an effective shear rate 



leff 



Viji + ?7272 



(10) 



Note that the second-order term does not depend on the 
elastic properties of the interface. 

The special feature of Eqs. (O and (fTU)) is that the 
effective shear rate j^ff felt by the interface differs from 
the applied shear rate 7 defined in Eq. ([1]) . Although the 
latter is set by the geometry of the system, the former is 
a dynamical quantity in the sense that it depends on the 
viscosities of both fiuids. However, jeff and 7 cannot be 
tuned independently since the continuity condition ([6]) 
for tangential forces requires 77171 = 77272 (see App. |B|). 
and then 



Li L2\ L 
— + — 7e# = - 7 



(11) 



This relation implies that the effective shear rate can 
in principle be tuned in the range < jgff < 2 7 by 
adjusting the experimental conditions. It is only when 
Vi — '72 that both shear rates coincide. 

Thermal fluctuations are accounted for through the 
white noise ip{q,t). This contribution stems from the 
random part of the stress tensor (see Appendix [C]) . We 
flnd that <^(q, t) simply follows the equilibrium distribu- 
tion with zero mean value and 



{^{<l,tMq',t')) 



ksT 
2fjq 



d{t-t'){2TT)^S{q^ 



[12) 
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Even though the system is driven far from equilibrium, 
there is no couphng between the noise and the external 
flow (at least up to ^(e^)) (H- 

At this point, it should be mentioned that corrections 
similar to Eq. ([9]) have been suggested in the context of 
sheared smectic phases (see [25| and references therein), 
or in the description of coarsening under shear [2^ . [2^ . 
The argument commonly invoked is the following. Since 
the deformation h{x, y, t) is dragged along the x direction 
by the external flow v*^"^ = w'^^^e^,, the time derivative 
involved in the equation of motion can simply be replaced 
by the total derivative 

dth-^dth + v^°^d^h , 

where v^^'^ is evaluated aX z = h. Of course if the slope 
of the velocity is continuous, i.e. if "qi = there is no 
ambiguity: the total derivative then reads dth + -yhdxh, 
yielding to the convolution integral in Fourier space. 
However, the situation is more subtle when rji ^ rj2- 
With the same argument, the velocity evaluated at z — h 
would be w*^"^ = 7i/i if /i > 0, whereas it would take 
the value = 72/1 if /i < 0. Clearly, there is an in- 
consistency in the reasoning. The only way to remove 
the indeterminacy is to follow the procedure detailed in 
Appendix |B] We emphasize in particular that the var- 
ious second-order contributions in Eq. (|5]) ~ ([8]) are all 
equally important to give the perfectly symmetric rela- 
tion Eq. 



IV. NONEQUILIBRIUM FLUCTUATIONS 

NESS properties of the interface are then extracted 
from the equation of motion ^ together with the noise 
correlations (HH). Unfortunately, this class of nonlinear 
stochastic equation cannot be solved exactly. We are thus 
forced to restrict the discussion to moderate shear rates. 
But this does not mean that should remain small, 
as might be expected from simple scaling arguments. In- 
deed, our study reveals that the parameter that actually 
governs the steady state properties of the system is the 
combination of the two dimensionless parameters 



(Til 



(13) 



Depending on experimental conditions, a can be small 
even if 7X2 is not p^ . The analysis presented in the 
following is valid in the range < a < 1. 

Eq. ^ is solved using a perturbation theory presented 
in Appendix IdI We first discuss the steady-state corre- 
lation function S'(q, 7) defined as 



hm {h{ci, t)h{q',t)) = 5(q, j){27:fd{q + q') 

r— >oo 



(14) 



It is given at equihbrium by S'(q, 0) = ksT/ {a{q-^ + 1^ ^)). 
Under shear, the spectrum is modified according to 

Siq,^) = S{q,O)[l~a^Iiqlc)cos''0g + Oia^)] , (15) 
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FIG. 2: Correction to the nonequilibrium fluctuation spec- 
trum AS'(q, 7) = a/^ksTll) X AS(q, 7), in the direction of 
1. In inset we show I{qlc). 



with 9q the angle between the direction of shear e^; and 
the wave vector q. The function I{qlc) depends only on 
the norm q ~ |q| of the wave vector. Explicitly, we get 



1 



d SCOS0, 



-1 



s|/(|x- 



s|) 



7(^) + /(s)+/(|x-s|) 



(16) 



with X — |x|, s — |s|, 9s the angle between x and s, and 
f{x) = (1 -|- x'^)/x. This integral cannot be evaluated 
analytically; the result of the numerical integration is 
presented in the inset of Fig. [5] [2^ . 

Eq. (|15p is the second main result of this paper. It 
shows that the coupling is maximum in the flow direc- 
tion {9q — 0), while the spectrum is not affected in the 
vorticity direction (^g = f )• As expected, the result is 
invariant by inversion symmetry q^ <-^- —qx or qy <-> —qy 
The correction AS'(q, 7) — [^(q, 7) — ^(q, 0)| is plotted 
in Fig. [5] in dimensionless units. It should be noticed 
that even though all the wave lengths are affected by the 
flow, the spectrum is mostly affected when q ~ The 
correction then vanishes like q~^ for larger values of 9, 
whereas it scales as AS'(q, 7) ^ g'* in the limit g — > 0. 

The mean square displacement of the interface is 
then obtained from the sum over all modes {h'^){j) = 
(27r)^^ J S{q, j)qdqd9q. Notice that the integral that de- 
fines the roughness at equilibrium is actually divergent. 
Regularization is achieved by means of a microscopic cut- 
off a so that (/i^)(0) = kBT/{2TTa) x \n{ya). In NESS, 
we find that the fluctuations are strongly reduced by the 
external flow 

{h')ij)^{h'){0)[l-Ka' + Oia^)] . (17) 

The correction is quadratic in the control parameter a. 
In this expression, is a universal constant in the sense 
that it depends neither on the properties of the fluids nor 
of the elastic constants of the interface [2^. Moreover, 
it is independent of the microscopic cut-off as soon as 
a/lc < 10^2. Numerically, we get K « 0.087. 
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FIG. 3: Dimensionless correlation function 2na /{kBT) x 
C(r, 7) for different values of the control parameter: a — 
(solid line), a = 0.5 (dot), a — 0.8 (dash-dot), and a = 1 
(dash-dot-dot), a) In the direction of shear {ip = 0), the 
correlation length decreases, b) In the vorticity direction 
(<p = 7r/2), the correlation length increases. 

Finally, we focus our attention on the NESS correlation 
function C(r,7) = {h{r,t)h{0,t)), which is the Fourier 
transform of the structure factor. As can be seen in 
Eq. (fT5|) . the spectrum S'(q, 7) is not modified for the 
modes that are perpendicular to the flow (i.e. wave vec- 
tors with = 0). However, performing the sum over all 
wave vectors gives rise, in direct space, to modifications 
even in the vorticity direction. Denoting ip the angle be- 
tween the direction of shear and the position vector 
r, the correlation function reads 

C(r, 7) = ^ [Ko (r/lc) - o? cos' ip {r/Q 

-a' sin' ^Cy (r/lc)] , (18) 

with Kq the modified Bessel function of the second kind. 
The corrections can be written Cx (x) — Ci {x) — C2 (x) jx 
and Cyix) = C^ix)^ where for i = 1, 2 

C.(x) = /°°ds4^i^X(s)s2-* . (19) 

Here, J„ is the Bessel function of the first kind. 

At equilibrium, C(r, 0) is proportional to fCo (tI^c) and 
behaves like C(r,0) oc exp[— r//c] at large separations 
r ^ If.. In other words the capillary length is also the 
equilibrium correlation length. In nonequilibrium situa- 



tion, the correlation function presents a number of inter- 
esting features. We first consider the direction of shear 
(y9 = 0. It clearly appears from Fig. [3la) that C(r, 7) 
decays faster that its equilibrium counterpart. This re- 
flects a decrease of the correlation length as the control 
parameter a is increased. The situation is different in 
the vorticity direction if — tt/2, as shown on Fig. [3lb). 
In this case the decay of C(r, 7) is slower than the equi- 
librium correlation function, indicating an increase of the 
correlation length with increasing shear rate. 

Since we do not have an explicit functional form of the 
correlation function, it is difficult to be more quantitative 
regarding the correlation length. The only conclusion 
that can be drawn from Fig.[3]is that the decay of C(r, 7) 
at large distances is not exponential anymore. Finally 
note that at intermediate angles {0 < ip < 7r/2), the 
correlation length can be either larger or smaller than Ic 
depending on the shear rate. 



V. DISCUSSION 

Complex fluids under shear is an challenging class of 
nonequilibrium systems. The simple system considered 
in this work is a fluid interface exposed to a Couette 
flow. First, it was shown that the time evolution of the 
deformation satisfles a mode-coupling equation. Note 
that equations similar to Eq. ([9]) have been suggested 
in the context of soft surfaces under shear (see for in- 
stance [lEHEIiSl)- a matter of fact, it belongs to the 
general class of the Kardar-Parisi-Zhang equation [s^l 
whose derivation is usually based on phenomenological 
grounds. Here however, Eq. ([9]) has been rigorously de- 
rived from hydrodynamics with no other assumption that 
inertial effects can be neglected. We now discuss our re- 
sults in view of recent studies. 



A. Comparison to experiments and simulations 

We have shown that interfacial fluctuations are 
strongly reduced by the flow. According to Eq. P7)l . the 
reduction is governed by a universal parameter: we pre- 
dict K « 0.087. Suppression of thermal capillary waves 
was indeed measured by Derks and collaborators in a 
recent experiment (l^ . [sij . This was achieved by using 
a phase-separated colloid-polymer mixture whose inter- 
face is characterized by a very low surface tension. The 
authors have used two compositions but the dynamical 
parameters are only available for the flrst sample (sample 
A, very close to the critical point). The thickness of the 
fluid layers are Li = 50 fim and L2 = 350 fim. Equilib- 
rium parameters obtained from the correlation functions 
are cr = 2.5 nN/m, Ic — 2.6 /im, and Tc = 13 s [l3 |. 
The viscosities extracted from the velocity profiles are 
rji — 7.8 mPa.s and 772 = 5.2 mPa.s [3l|. Unfortunately, 
there are only two experimental points (out of 6) that 
correspond to the condition a < 1. To fit the data for 
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FIG. 4: Mean square displacement of the interface a a func- 
tion of the control parameter a. The solid circles refer to the 
experiment of Ref. [ij] (sample A). The dotted and solid lines 
are the theoretical predictions of Eqs. (|17|l and (|20|l , respec- 
tively, with Ktheo ~ 0.087. The dashed line corresponds to 
the fit of Eq. with Kfit ^ 0.246. The dash-dotted line is 
the model of Derks et al. 14] . 



higher values of the control parameter, we assume that 
Eq. (fTT)) can be extended according to 



(ft2)(7)^(/i^)(0) 



1 



1 + Ka^ 



(20) 



We then perform a linear regression of the quantity 
(/i2)(0)/(/i2)(7) as a function of a^. We find Kfit « 0.246, 
almost 3 times the expected value. The comparison be- 
tween our theoretical predictions and the experimental 
measurements is shown on Fig. S) As can be seen the 
agreement is only qualitative, even though it is difhcult 
to be really conclusive since most data are outside the 
vahdity range of Eq. ((T71) (3^ . 

In Ref. |14i | , the suppression of thermal capillary waves 
is supported by a remarkably simple model. The authors 
state that only the slow modes with '^Tq > 1 are affected 
by the flow. Solving this inequality defines to two wave 
vectors qi and q2. It is then assumed that wave vectors 
with qi < q < q2 no longer contribute to the fluctuation 
spectrum, whereas the contribution from wave vectors 
outside this range is left unchanged. Their prediction 
of the interfacial roughness as a function of the applied 
shear rate 7 is in very good agreement with experimen- 
tal data - see Fig. [H Still, we argue that the hydro- 
dynamic analysis presented in this paper provides new 
insights into the problem. In particular, we have shown 
that fluctuations are driven out of equilibrium by the ef- 
fective shear rate 'y^ff rather than the applied shear rate 
7. For sample A presented above, the ratio of effective 
to applied shear rate is jeff/j ~ 0.8. Now, suppose that 
the experiment is done with the same sample composi- 
tion excepted that the thicknesses of the fluid layers are 
inverted to Li = 350 /.tm and L2 — 50 fiTa. In this case, 
the ratio would be 'jeff/j ~ 1.2: for the same value of 
the applied shear rate, a is thus expected to increase by 
50 %. The ensuing difference regarding the mean square 



displacement should be substantial enough for the dis- 
tinction between applied and effective shear rate to be 
experimentally observed. This would provide a direct 
validation of our analysis. 

Finally, let us discuss the correlation functions. Our 
study reveals that the static structure factor is strongly 
modified in the flow direction, whereas fluctuation modes 
in the vorticity direction are not affected. This result 
confirms similar conclusions obtained from molecular dy- 
namic simulations [33|- Coming back in real space, we 
find that the correlation length decreases in the direction 
of shear. This conclusion agrees with Monte Carlo sim- 
ulations of a driven interface [lB|- Physically, it would 
mean that the external flow acts as an effective potential 
whose strength should increase with the shear rate. How- 
ever, these findings are in disagreement with the results 
of Derks et al. that observe an increase of the correlation 
length [T5| • Note that experimental data are fitted using 
the equilibrium correlation function. This procedure is 
mainly based on the short distance behavior of C(r, 7). 
In contrast, our conclusions are drawn from the asymp- 
totic behaviour of the correlation function. It is thus 
possible that the discrepancy comes from how the cor- 
relation length is defined in practice. Other possibilities 
are discussed below. 



B. Conjectures and future directions 

The equation of motion ([§]) is a direct generalization 
of the dispersion relation for capillary waves. In par- 
ticular, it is based on the same assumptions: continuum 
hydrodynamics and small-gradient expansion. The perti- 
nent question is whether the apparent mismatch between 
theoretical predictions and experimental data can be the 
consequence of another mechanism [3^ . In fact, several 
explanations are plausible. 

Firstly, we remark that the diameter of the colloids is 
of the order of 140 nm. The capillary length is thus only 
ten or twenty times larger than the "microscopic" cut- 
off, and it is not clear whether the separation of length 
scales assumed to get the universal quantities {K, X, . . .) 
is really achieved for this sample. 

Secondly, Monte Carlo simulations of Smith et al. re- 
veals that the average width of the interface is signifi- 
cantly reduced by the flow [l^. In the capillary wave 
theory, it is assumed that the density proflle behaves like 
a step function with a strictly vanishing width. But it 
is known that physical parameters are directly related to 
the shape of the true density profile [33^. For instance, 
the surface tension is given by 



dz(|^ 

oz 



(21) 



with p(z, 0) the equilibrium profile and m a constant pro- 
portional to the second moment of the intermolecular po- 
tential [sl] . In the small but finite interfacial region, the 
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density is varying rapidly so that the fluids are highly 
compressible. The density should therefore be modified 
by the flow and become a function of the applied shear 
rate p(z,j). A plausible consequence is that the bare 
surface tension might depend on 7 as well. According 
to Eq. (|2ip. the thinning of the average profile observed 
in hj is expected to give rise to an increase of the sur- 
face tension, and then to an additional reduction of the 
fluctuations. 

The latter assertion is certainly hypothetical but we 
argue that the measurement of the correlation function 
in the vorticity direction could be conclusive. Indeed, an 
increase of the bare surface tension would lead to an in- 
crease of Ic and thus of the NESS correlation length. But 
a modification of C(r, 7) arising from a{j) is expected to 
be isotropic, whereas the hydrodynamic theory predicts 
a modification that is anisotropic. This question clearly 
deserves to be reconsidered both from the experimental 
and theoretical viewpoint. 
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APPENDIX A: DERIVATION OF THE 
BOUNDARY CONDITIONS 

Since the analysis presented in this work involves non- 
linear contributions, particular care has to be given to 
the conditions at the interface. In this appendix, we first 
remind the derivation of the kinematic relation Eq. ([5]). 
We then focus on the stress continuity condition Eq. ([6]). 



1. Kinematic condition 



VI. CONCLUSION 

Let us briefly summarize the main results of this paper. 

(i) We have generalized the hydrodynamic derivation of 
the dispersion relation to include the effect of an applied 
shear rate on the dynamics of capillary waves. Eq. ^ 
reveals that the system is driven out of equilibrium by an 
effective shear rate that differs from the applied one. The 
mode-coupling structure of the equation also shows that 
nonequilibrium fluctuations cannot be described in terms 
of a simple rcnormalization of the interfacial tension. 

(ii) The analysis of the mode-coupling equation shows 
that the fluctuations arc smoothed out by the flow, in 
qualitative agreement with experiments and simulations. 
We find that NESS fluctuations can be described in terms 
of a single control parameter defined in Eq. (|13p . Conclu- 
sions regarding the height correlation function are more 
questionable. If the theory and the simulations predict 
a decrease of the correlation length in the direction of 
shear, the reverse is actually observed on the experimen- 
tal side. The discrepancies might be due to the finite 
thickness of the interfacial region and its response to the 
external fiow. Molecular dynamics simulations of inter- 
faces under shear should help to clarify this point in a 
near future. 

In conclusion, we have derived a hydrodynamic theory 
to account for the coupling between a shear flow and the 
fluctuations of an interface. Previous investigations on 
driven interfaces have revealed partial agreements but 
also signiflcant discrepancies between experiments and 
numerical simulations. Although our work brings new 
insights into the issue, it seems that further investiga- 
tions will be necessary to definitely reconcile experiment, 
theory and simulations. Ultimately, achievement of this 
program should allow to reach a deeper understanding of 
soft matter systems under shear. 



It is convenient to represent the interface with the func- 
tional f{x,y,z,t) = z — h{x,y,t). The unit vector n nor- 
mal to the surface is given by 



n = 



V/ 

|V/| 



(Al) 



Let us focus on a point M(x, y, z) that belongs to the 
interface. It satisfies f{x,y,z,t) — and then 

d/ = dtfdt + V/ • dr = . 

If w is the velocity of the point M, the infinitesimal dis- 
placement reads dr — wdt and therefore V / ■ w = —dtf- 
Together with Eq. (|XT|) . we then obtain 



|V/| 



Identifying the normal velocity of the interface w • n with 
the normal velocity of the fiuid v • n, one finally gets 



dth 



V • n 



(A2) 



Here, we have also used the fact that dtf — —dth. This 
relation together with the definition ([7]) directly leads to 
Eq. ^. 



2. Dynamical condition 

To derive the stress condition ([S]), we consider an in- 
finitesimal surface element S of the interface, bound by a 
closed contour C. The contour is chosen circular with ra- 
dius e. We also define a small cylinder V, centered on the 
contour C and with height 2e. The forces acting on the 
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volume element V are the body force density f , the ten- 
sion force exerted along the perimeter C, and the surface 
force exerted by the fluids above and below the cylinder. 
We discuss separately each contribution. 



a. Tension force 

We first discuss the restoring force exerted on the con- 
tour line C. We assume that the contour is travelled 
counterclockwise. The (local) orthonormal basis associ- 
ated with the line is denoted {m, s, n}, with m the tan- 
gent vector, s the vector normal to C but tangent to 5, 
and n the vector normal to S pointing towards the upper 
fluid. In this representation, the external force exerted 
on the contour C reads 

Ft = — <j) (Tsdl = ("^ X n) dl ~ a (j) {dl X n) , 

with dl = mdL Then, according to the Stokes relation 
^ dl X A = J (dS X V) X A, one gets 



Ft = CT y d5 (n X V) X n , 



with the surface element dS = dSn. 

Finally, using the general identity (A x V) x A = 
i V (A • A) - A (V • A) and because n • n = 1, we find 



Ft^-a J dS (V • n) n 



(A3) 



b. Surface force 

The component r^,y of the stress tensor T corresponds 
the /I*'' component of the force (per unit area) acting on a 
surface element normal to the direction v. In particular, 
the force density exerted by a fluid on a surface element 
with outward normal n is T • n. 

Let us apply this definition to the elementary cylinder 
of volume V. The cylinder is sitting across the surface, 
the extension being e on each side. Regarding the upper 
surface of the cylinder, the stress exerted by fluid 1 is 
characterized by the (outside) normal vector rii = n so 
that the force density reads Ti • n. Similarly on the lower 
surface one has n2 = — n and the stress is — T2 • n. The 
total surface force is then 



JdSiTi 



T2) 



(A4) 



Now if e is the typical lengthscale of the volume element 
V, then the surface contributions scale as whereas the 
acceleration and the body forces scale as e^. Hence in 
the limit e — > the latter can be neglected and surface 
forces must balance. However, since the surface element 
is arbitrary, the integrand must vanish identically so that 
one flnally gets the desired result 



(Ti-T2)-n = cr(V-n) n 



(A5) 



Finally notice that the divergence of the normal vector 
reads in the Monge representation [20| 



V ■n = -V^h + 0{h^) 



(A6) 



APPENDIX B: NONLINEAR RELAXATION 
EQUATION 

The aim of this section is to derive the equation of 
motion ([9]). The fluctuation modes of the interface are 
expected to be deformed by the shear flow. As we shall 
see, the description of this effect requires to go beyond 
the usual linear analysis. We focus here on the relaxation 
dynamics in the absence of noise. Thermal fluctuations 
are treated in Appendix [Cl 

To begin with, we recall the Stokes equation 



?7,VV - Vj3 + Pig 



(Bl) 



with g = —gSz the gravitational acceleration. In the 
small-gradient approximation, we can express the defor- 
mation profile as h(x,y,t) = su{x,y,t) with u{x,y,t) ~ 
0(1). The dimensionless parameter that governs the 
small-gradient expansion is 



all 



(B2) 



We then assume that the fields that satisfy the Stokes 
equation can be expanded in powers of e 

F{z)=F'^''\z)+eF'^^\z)+e'F'^''\z) + 0{e^) , (B3) 

where F stands for v, t or p. (For clarity reason, only 
the z-dependence of the fields is kept explicitly in this 
appendix) . The solution at order e° satisfies Eq. (jBll) ; it 
corresponds to the Couette flow solution for a flat inter- 
face. The other terms of the expansion (|B3[) are solution 
of the following equation 



Vp("' = , Vn > 1 



(B4) 



c. Force balance equation 

The force balance on the volume element V enclosing 
the surface S deflned by the contour C then reads 



Ft 



dV^f . 



At order e , the solution in the absence of shear would 
lead to the usual dispersion relation for capillary waves. 
Here, the calculations are carried out up to order e^. 

Due to the linearity of Eq. (|B4p , the coupling between 
different orders only occurs through the boundary con- 
ditions at the interface. This suggests to use a recursive 
method in order to solve the problem [23l| . 
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1. Solution at order flat interface 

For a flat interface, the solution is the well-known so- 
lution for the Couette flow: 



(B5) 
(B6) 



with i = 1 if z > and i — 2 if z < 0. is a constant. 
The stress tensor then reads 



77,7, 
t(")(z) = ( 



(B7) 



In particular, the continuity relation ^ for the tangen- 
tial stress implies r/171 = 77272- 

Note that we have chosen Li, Vi, L2 and V2 so that 
the plane of zero shear is the horizontal plane z = 0. At 
constant shear rate, any other situation can be deduced 
thanks to a simple Galilean transformation. 



2. Solution at order e^: linear relaxation 



The first-order solution satisfies Eq. (|B4p with bound- 
ary conditions that involve the solution at order e°. Ac- 
tually, there are two kinds of conditions that needs to 
be enforced: continuity conditions at the interface and 
limit conditions far from the interface. Strictly speaking, 
the latter have to be expressed at the boundaries of the 
shear cell, z —> Li or z ^ where Li and L2 are 

macroscopic lengths. In this work, we are interested on 
thermal fluctuations that develop on microscopic length- 
scales (i.e., with wavelength much smaller than Li or L2). 
We can therefore assume that for any field F, the terms 
of the series (|B3p with n>l vanish at infinity 



(B8) 



lim = , Vti > 1 

z — '±00 



In the following, we first derive the boundary conditions 
at the interface and then solve the problem in Fourier 
representation. 



a. Boundary conditions 

To express the continuity conditions, let us define for 
a field F the notation 



[F]^=F{zt)-F{z^) 



(B9) 



with F{z^) = limz-tzf, F{z), the limit being taken for 
z ^ zq. The first condition that needs to be enforced is 
the continuity of the velocity at the interface 



Using the expansion (jB3p together with the Taylor ex- 
pansion of w{h — eu), we obtain at first order 







From the result (jB5|) we get 9^v('^^ = jtex, and finally 



,(1) 
.(1) 



= (72 -ii)u , 







= 



(BlOa) 
(BlOb) 
(BlOc) 



Next, we have to enforce the condition of stress conti- 
nuity Eq. At first order, we obtain 



with fj, — x,y or z. Following the result (|B7p for the 
stress tensor, we find 

dzt^^Kz) = Pzg^i^z ■ 

The continuity conditions for the stress tensor then reads 
= , (Blla) 



''xz 



i-yz 



''zz 



= 



= Apgu — a\7 u , 



(Bllb) 
(Bllc) 



with Ap — P2 ~ Pi- 



b. Method and solution 

The solution of the problem is expressed using the two- 
dimensional Fourier transform 

F{q,z,t) — J d^r cxp[— iq • r]i^(r, z, t) , 
F(r, z, i) = y exp[-iq • r]F(q, z, t) , 



(B12) 



with r = {x,y) and q = {qx,qy)- The difficulty of the 
algebra comes from the fact that we need to solve the 
full 3-dimensional problem. It then appears judicious to 
define a new orthogonal coordinate system that would 
account for the symmetries of the problem. To this aim, 
the vector fields are decomposed into their longitudinal, 
transverse and normal components [sgI [37j . This defines 
a new set of orthogonal unit vectors {l,t,ez) where 1 is 
the unit vector parallel to q, and t the in-plane vector 
perpendicular to 1 and e^. These vectors are expressed 
in the cartesian basis {ex,ey,ez) as 



1 



qx 



1y. 



qy 
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The velocity is then written as v^^^ = v 
Inserting this expression in Eq. (|B4p leads to a system of 
differential equations for the Fourier-transformed quan- 
tities 



together with the divergenceless condition 
iqvl^^ + 9^4^) = . 



(B13a) 
(B13b) 
(B13c) 

(B14) 



It is not difficult to show that the solution of Eqs. (|B13|) 
and (|B14p can be written 



vPiz) = A^^'e-"^^'^ + Bl'>qze 



l(i) 



-<l\^\ 



,(1) 



(z) = ^e-«l^l , 



(B15a) 
(B15b) 
(B15c) 



Ap\ and c|^^ are the (yet undetermined) integra- 

tion constants, the subscript i taking the value i = 1 if 
z > and i = 2 if z < 0. (There should be a priori 
no confusion between the subscript and the imaginary 
unit i). The solution for the longitudinal component fol- 
lows directly from (|B14|) . Note that the condition ([BS]) 
has already been taken into account in (jBlsp . 

Next, we have to express the boundary conditions 
Eqs. (|B10[) and (jBll[) in Fourier space and in the new 
system of coordinates. The condition (jBlOcp for the nor- 
mal component is straightforward 



= , 



(B16) 



where of course vi^\q,z,t) now stands for the Fourier 
transformed quantity. Then projecting conditions (|B10a[) 
and (jBlObp for the parallel components onto the direc- 
tions 1 and t respectively gives 



and 



„(i) 



.,(1) 



-(72 



-(71 



71 ju 



■72)" 



(B17) 



(B18) 



Regarding the continuity of the stress, the condi- 
tion (|Bllc[) is now expressed as 



-I- 2r,,a,«(i) - (g2 + u , (B19) 
Jo 



where we define the capillary length by — a/{Apg). 
Finally, the projection of (jBllap and (jBllbp leads to 



iqv 



(1) 



= 



(B20) 



and 



(1) 



= 



(B21) 



With the boundary conditions ((BT6)) - (|B2l|) . we can 
now evaluate the 6 integration constants. We find 



(1) 



A. 



(1) 



Ariq 



a (q^ + 



B 



(1) ^ ^(1) 



(1) 



B. 



2riq 

4(1) , ■ 1^ ■ f \ 
+ 177^12 im -V2)U 
2r]q 



2'qq 



7/1 ) U 



C. 



(1) _ 1y 



2r]q 



(B22a) 
(B22b) 
(B22c) 
(B22d) 
(B22e) 



with f] — + 772)72. 

Although the first-order solution already depends on 
the shear flow, it is claimed in the main body of the paper 
that the calculations have to be performed up to second 
order. To understand this point, consider the kinematic 
condition ([5]). Since v'^*'^(0) = 0, we can write at first 
order 



a^/i = (0) 



1 



-h , 



with Tq = —I/A'^K It appears from (jB22aP that the 
shear rate is not involved in the equation of motion at 
linear order. The calculations have thus to be extended 
to include the first nonlinear contribution. 



3. Solution at order e^: mode-coupling equation 

a. General solution 

The solution at second order satisfies the same set of 
linear differential equations (|B13p and (|B14p . so that it 
takes the form 



(z) = Af)e-''l^l+sf gze-'l^l 
.r^(z) = cf)e-N, 



(B23a) 
(B23b) 
(B23c) 



Again, is deduced from the incompressibility condi- 
tion (|BT4| . 

It appears that not all the integration constants are 
needed for our purpose. Indeed, let us express the kine- 
matic condition ^ up to second order. In direct space, 
it reads 

edtu + e^ud^ud.v^^^ (0) + Vp • v[|^^ (0) 



(0) + e\d,vi^\0) + e\f\0) 
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where we have used the fact that v(°)(0) = 0. Note that 
the passage from ^ to this latter equation leads to an in- 
determinacy. Indeed, although the velocity field v is con- 
tinuous at z = ft., this is not the case for each term of the 
expansion (|B3p at z = 0. In this appendix, we thus as- 
sume that the limit z ^ is taken from above. It can be 
checked that the same results are obtained if one would 
take the limit from below. We now switch to Fourier 
space, where the product of two functions f{r).g{r) is 
transformed into the convolution product 



if* 9) (q) 



d^k 



/(k).9(q-k) 



(B24) 



(2^)2 

We then have in Fourier representation 

edtu — ev'^J'^ (0) + [u* dzV^p (0)) — 171 [u * q^v^ 

-i(vW(0)*g,u) -i(i;W(0)*gj,u)] +e''vf\Q) . 

With the results of the preceding section, the equation 
of motion can now be written as 



edtu — eu — ie'^ji (u * Qxuj 



(B25) 

(2) 

Therefore, only the integration constant Al is required 
to get the equation of motion. 
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b. Integration constants 

At second order, the continuity condition ([5]) for the 
z-component of the velocity is expressed in direct space 
as 

,2 



leading to, in reciprocal space. 











,(2) 



i (72 - 71) (" * Qxu) 



(B26) 



(2) (2) 

This relation actually involves both A\ and A2 ■ To ob- 
tain one last equation, we still need to express the normal 
stress condition 



[n-t-n]^ 



Together with (|B7p and (jBlip . it reads at order 



= 



This result is finally expressed in Fourier representation. 
The derivative of the stress tensor at first order is then 
obtained from Eq. (|B13cP 



where we have used (IB14P to get the last equality. From 
the continuity relation (|B20p . it is not difficult to obtain 
the condition 



(2) 







that leads to the simple relation 



m^?^ + V2A 



(2) 







Inserting this result in Eq. (|B26p . we finally get 



^'1^ = * q^u) 



A. 



(2) 



i"/2[u * q^u) 



2fj 

m - vi 

2f] 



(B27) 

(B28) 

(B29a) 
(B29b) 



If we factorize all second-order contributions in (|B25p . 
we then arrive at the equation of motion 



dtsu 



1 . 2?7i7i / s 
— eu — I [eu * qx£u) 



(B30) 



Coming back to the deformation field h = eu and using 
the fact that 77171 = 77272, we can write the equation of 
motion in its definitive form Eq. 



APPENDIX C: FLUCTUATING 
HYDRODYNAMICS 

In this appendix, we consider the Landau-Lifshitz 
equation of linear fluctuating hydrodynamics [l3|. We 
thus need to adapt the perturbative scheme developed 
in Appendix [B] to the Stokes equation ^ including the 
random part of the stress tensor s. To this aim, it should 
be noticed that the noise correlations (s^iyS^^v) given in 
Eq. ^ scale as fc^T and are therefore proportional to e"^ 
[see Eq. (IB2p ]. From this viewpoint, the random part of 
the stress tensor can be identified as a first- order contri- 
bution. 

We then follow the same recursive method as in Ap- 
pendix [B] There is no modification at order e° and we 
now discuss the first and second orders of the expan- 
sion (lB3l). 



1. Solution at order 

At first order, the stochastic version of the Stokes equa- 
tion ([2]) reads, in the (l,t,ez) basis, 

77, {dl - q^) vf^ = 7gp(i) - iqsii - d^su , (Cla) 

(Clb) 



rj,{d^^-q^)vl^>^-d^s,^, 
77, - 92) yW ^ _ ,qs,^ _ 0,s,, , (Clc) 

together with the incompressibility condition 

iqvl^^ + d.v^p = . (C2) 



12 



As explained in the text, the fluctuations of the random 
forces in the bulk are not affected by the shear flow (at 
least in the regime considered in this paper). The corre- 
lations are then given by the fluctuation-dissipation the- 
orem [l3, [3 ■ In this representation, we have 

X {2Tif6 (q + q') 5{z - z')S{t ~ t') , (C3) 

where the subscripts stand for I, t or z. 

The solution of Eqs. (jCip and (|C2p is then obtained 
using the following Green's function identity [l3| 



1 

2,' 



'1 = -6{z-z') 



Regarding the boundary conditions, the stress tensor t 
in (|B10[) and (jBlip has now to be replaced by the total 
stress tensor T = t + s [l^, [l^ . After some algebra, we 
obtain the evolution equation for the interface 



dMq_, t) = - — /i(q, t) + (^(q, t) 



(C4) 



where the noise (/3(q, t) is given by 
4?7 Jo 



+ 



4:1] J_ 



- sii + 2isiz) 
sii - 2isiz) . 



(C5) 



Together with the correlations (|C3p . we finally get the 
result announced in Eq. p^ . Note that Eq. (|C5p gen- 
eralizes the result of Grant and Desai for a liquid-gas 
interface (in the absence of shear) [l^ . 



2. Solution at order e 

At second order, the velocity field satisfies the following 
set of equations 



together with the continuity equation 



, 







(C6a) 
(C6b) 
(C6c) 



(C7) 



The noise enters the problem only through the bound- 
ary conditions and the first order contribution. Again, 
the boundary conditions are the same as in Appendix [Cl 
excepted that the stress tensor t is replaced by the total 
stress tensor T = t -I- s. 

The calculations are pretty lengthy so that we only give 
the final results. The equation of motion is obtained from 
the kinematic condition Eq. ([5]). Up to second order, we 
find 



dth = h - ijeff [h * {qxh)] +ip-q - [(q^i) *h] - (qx [ (q02) * h]j • , 



(C8) 



r 



where 4>i and (/)2 are defined as 
1 f^°° 

<j)i{q, t)^ — dz (1 - qz) e «^ (s^^ - su -f 2isi,) 

4m Jo 
1 /-o 

+ dz (1 -I- qz) e"?^ {s,, ~ su - 2isu) , 

(C9) 

and 

1 /■"^°° 

(j)2{q,t) = — / dze {sit + istz) 
Jo 

4--^ / dze'^^sit - tstz) . (CIO) 
2??'? J -00 

It can be noticed that the additional terms involving 
and 02 are not induced by the external flow. At equi- 
librium, these second-order contributions would actually 
renormalize the relaxation dynamics of the interface. In 



nonequilibrium situation, it is argued in the next ap- 
pendix that those terms do not contribute to the fluc- 
tuation spectrum at the order considered in this work. 
Therefore, we do not to mention them in the main docu- 
ment. This choice is made in order to improve the read- 
ability and to focus the discussion on the main outcomes. 



APPENDIX D: SOLUTION OF THE 
MODE-COUPLING EQUATION 



1. Linear noise 

Eq. ([9|) is solved using a perturbation theory. To this 
aim, let us define the time Fourier transform /(q, cj) of 
any function /(q, i) as /(q, w) = / dtexp[—iujt]f{q,t). 
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The mode-coupling equation can then be rewritten 

X / dkk^h{k,n)h{q-k,uj , (Dl) 



with dk = d^kdri/(27r)'^. The bare propagator R{q,uj) 
is given by 



1 + iUTn 



(D2) 



The solution of the problem is then obtained as the solu- 
tion of a Dyson equation. The calculations are performed 
up to second order in the small parameter a (which 
is proportional to jeff)- Using the shortland notation 
/iq = /i(q, w), we get 



Rqipq- ijeffRq / dkA;j;i?k-Rq-k<^k<^q-k " 7ejff-Rq J dkdk'/Cj; fc^i?ki?q-k^k' -Rq-k' <^q-k<^k"^q-k' 

- -Tgjjr^q / dkdk'/Cj;fc^^k^q-k^k'-Rq-k-k'',5k'^k''^q-k-k' + O , (D3) 



From this expansion, the evaluation of the correlation 
function {hqhq>) is pretty lengthy but presents no con- 
ceptual difficulty since the four-point correlation func- 
tions of the noise are given by the Wick's theorem. The 
result is then Fourier-transformed back in time represen- 
tation, leading to Eqs. (fTS]) and (fT6|) . 

2. Nonlinear noise 

The same scheme is used to treat the nonlinear noise 
contributions. First, we need to specify the statistical 
properties of (pi and (p2- From (|C3p . it is not difhcult to 

get {(f>i) = {(1)2) ~ 0, and for v ~ 1 or 2, 

(D4) 

and 

(0^(q,i)¥>(q',t')) = O. (D5) 



The solution is then expressed as a Dyson equation sim- 
ilar to Eq. (|D3p . Although the expression of hq involves 
more terms, it can be shown that most of the additional 
contributions to (hqhqi) actually vanish (either because 
they involve an odd number of noise terms, or because of 
the vanishing cross-correlations between ip, (pi and 02)- 
As a matter of fact, additional shear-induced contribu- 
tions are at least of order , and thus do not contribute 
to the modification of the spectrum at the order consid- 
ered in this work. 

Still, the nonlinear noise does give rise to non- vanishing 
contributions at order e*. But those terms are equilib- 
rium contributions since they do not involve the shear 
rate. Consider for instance the mean square displace- 
ment: those terms would simply renormalize the equilib- 
rium value (/i^)(0). This is why we have carefully plotted 
the ratio {h"^) {j) / {h'^) (0) , that represents the modifica- 
tion of the fluctuations by the shear with respect to the 
equilibrium situation [38| . 
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